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J>^' Abstract 

| We propose extended version of stationary string model to describe the nuclear 

attenuation. This model takes into account flavour content of particles and allows 
to include into consideration all hadrons created from string. The predictions of the 
model are compared with experimental results obtained by HERMES collaboration 

' (DESY) on different nuclei (N and Kr). 
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<N ■ 1. Introduction 
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C\l ■ The study of the hadrons production in deep inelastic lepton-nucleus scattering 

offers the possibility to investigate quark propagation in dense nuclear matter and the 
space-time evolution of hadronization. In particular, the measurements of high energy 
hadrons attenuation in nuclear matter is a well-known tool to specify the parameters of 
models related to the early stage of particle production. 
^ . Hadron production in deep inelastic scattering (DIS) of leptons on nucleons goes through 

>• ! two stages. On the first stage the virtual photon knocks out the (anti)quark from nu- 

cleon, which on the second stage produce a observable hadron. The strong interactions 
theory - perturbative Quantum Chromo Dynamics (QCD) at present can not decribe the 
process of quark hadronization because of in this case "soft" interactions play a greater 
role. Therefore, experimental investigation of quark hadronization and building of phe- 
nomenological models describing the most common characteristics of given process is of 
importance for development of QCD. 

Let us consider the simplest possible state, a color-singlet quark-antiquark system. The 
lattice QCD studies lend support to a linear confinement picture, i.e. the energy stored 
in the color dipole field between a charge and an anticharge increases linearly with the 
separation between the charges. This is quite different from the behaviour in QED and 
is related to the presence of a triple-gluon vertex in QCD. If the tube is assumed to be 
uniform along its length, this automatically leads to a confinement picture with a linearly 
rising potential. From hadron spectroscopy the string tension, i.e. the amount of energy 
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per unit length is deduced to be k ~ 1 GeV/fm |T7|]. The linear confinement provide a 
simple explanation for the existence of linearly rising Regge trajectories. In this scheme 
the string tension determined by the Regge trajectory slope: k = (27ra R ) , which is also 
~ 1 GeV/fm Eg,E§. 



Hadronization process takes place at distances of few fm, right after the deep inelastic 
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interaction with nucleon. As this corresponds to the size of the nucleus, it is clear that 
answers to these questions can be obtained through investigating the processes of hadron 
production in DIS of leptons on nuclear targets. The available experimental results were 
obtained on different nuclear targets with (anti)neutrino [[[], |2|, muon H |], electron |J 
and positron || beams. The experiments showed the presence of nuclear attenuation 
effect for charged hadrons with xf > 0.1 (xf is Feynman's variable), which is essential 
in the energy range of ~ 10 GeV and almost vanish at energies of more than 50 GeV. 
It means that the multiplicity of hadrons on nuclei is different from the multiplicity on 
deuterium (both per nucleon). The experiments showed also the strong dependence of 
this effect on the energy transfered by leptons to nucleus, for electroproduction processes 
it is energy of virtual photon. 

In the meantime, experimental investigations let to development of phenomenological 
models. For description of nuclear attenuation, mainly the stationary string model (SSM) 
was used |10fl . It is supposed that after the DIS of lepton on internuclear nucleon, a 
color string is stretched between the knocked out quark and the nucleon remnant, a string 
which consists of gluons with constant tension over it's length. The color field of the string 
creates quark-antiquark and diquark-antidiquark pairs, which lead to breaking of original 
string into many short strings. This process results in creation of many string-hadrons. 
The quark which is on the fast end of string, while passing through the nucleus, can 
exchange the color with one of the quarks of the nucleon, lying within it's trajectory. As 
a result, the string breaks on two strings, the energy of the main string reduces and as 
a consequence, less fast hadrons will be produced in nucleus (per nucleon), than on free 
nucleon. It is also presumed, that energy and momentum of interacting quarks change 
insignificantly (color interaction). Some authors suppose that the color interaction cross- 
section is constant o c q = const || [TQ[| , others believe that it depends on Q 2 fl2||(Q 2 =-g 2 , 



where q - is the 4- momentum of the virtual photon). We return to this question later. 
The model supposes that the fast end of the string always contains a quark. It means 
that the validity of the model is limited by the valence region xsj > 0.2 (xBj is Bjorken's 
variable), because in the "sea" region, it is likely that an antiquark can be on the fast 
end of the string, and model does not describe the antiquark interaction. The quark may 
also deviate from it's primary movement direction as a result of interaction with nucleon 
by means of Pomeron or other Reggeon exchange. 

Let us also briefly consider other models for nuclear attenuation. In |T3j authors pro- 
pose a model of hadronization of highly virtual quark in nuclear environment by means 
of gluon bremsstrahlung and the deceleration of the quark as a result of radiative energy 
loss. This model practically transforms into the SSM at Q 2 < 3 GeV 2 . Recently in |14j] 
the authors supposed that in deeply inelastic eA collisions in the framework of multiple 
parton scattering, the quark fragmentation functions are modified due to higher-twist ef- 
fects. The hadronization of the quark takes place outside of the nucleus, which is correct 
in very high energy ranges. The model has lower limit of application over Q 2 , which is 
not clearly stated by the authors. Also, unrealistic nuclear density distribution functions 
for middle and heavy nuclei is used. These questions are important in sense of application 
to the HERMES kinematics. 

Conclusion is that in virtual photon energy region v = 5 - 25 GeV and <Q 2 > equal 2 
- 3 GeV 2 which are interesting for us in connection with the HERMES experiment ||, 
most suitable phenomenological model which can be applied for quantitative calculations 
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is SSM. Therefore, let us outline those points of this model which are improved in present 
paper: 

- description in final state only hadrons summed over kinds and charges; 

- the simple description of nucleon structure without " sea" partons is used. It leads to 
the limitation of applicability of the SSM only for valence region (xBj>0.2); 

- using only part of the quark-nucleon cross section connected with color interaction; 

- the problem connected with description of the Q 2 dependence of quark-nucleon and 
" hadron" -nucleon cross-sections, where "hadron" means colorless quark- ant iquark system 
on early stage of it creation without "sea". 

We propose development of SSM, including in consideration the kind and charge of 
hadron in final state. For this goal we take into account flavours and flavour contents 
of (anti) quarks and diquark from initial nucleon and quark- ant iquark pairs produced in 
color field of string and also of final hadron. 

Second innovation is that in framework of some fragmentation scheme we calculate the 
partial energy Zh = Ehjv (where Eh and v are energy of hadron and virtual photon in 
target rest frame, respectively) and constituent formation length for any hadron produced 
from string. It allows to calculate the nuclear attenuation with high accuracy. To simplify 
the calculations we limited ourselves only a three fastest hadrons. In other hand because 
of all calculations in this work were performed with the cut on Zh > 0.1-0.2 , then one 
can take in account the limited number of fast hadrons in the string. This is reasonable 
because of according to multi-periferical model the parton energy during hadronization is 
divided between the final hadrons in the following proportion: Zh\ ~ 1/2; Zh2 ~ 1/4; 2^3 
~ 1/8 and etc. 

We considered in final state only charged mesons, because the mechanism of (anti)protons 
has some peculiarities, consequently our model in present time does not claim to describe 
the nuclear attenuation of (anti)protons or h + and h~. We hope to include them in con- 
sideration in near future. 

1. Description of developed approach 

The proposed approach is based on the SSM, in which it is supposed that into the 
nucleus at the point with longitudinal coordinate z and the impact vector b the DIS takes 
place on one of nucleons (proton or neutron). Between the knocked out (anti)quark and 
nucleon remnant the color string is stretched . The maximal length of string is L~z//k 
(k is the string tension). Then the breaking of string takes place by means of quark- 
antiquark and diquark-antidiquark pair production in color field. The first constituent 
(anti) quark of the hadron with definite kind, charge and partial energy Zh is produced on 
a distance l c from the point of deep inelastic interaction. The l c is satisfied the condition 
< l c < lemaxi where l cm ax is the maximal distance on which the first constituent of the 
hadron may create, l cma x = (1 - Zhjvj'k .This relation was derived in || based on more 
general reasons than the string model . However this relation could be explained very 
simple in the framework of the string model if we suppose that hadron in own rest frame 
is the string with l = m^/k length, and in the system where the hadron energy is equal 
Eh, \=loEh/rrih and l c = L - 1, thus the relation is realized. 

Usually to describe the nuclear attenuation effect the normalized (per nucleon) ratio of 
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hadrons multiplicity produced on nuclear and deuterium targets is used. This ratio could 
be expressed via the ratio of corresponding fragmentation functions: 



. 2D\{z h ) 

where D^(^) and D\(z h ) are the quark fragmentation functions in deuterium (as- 
sumed same as in vacuum) and in A nucleus respectively. 



p poo pl c max " 

D h A (z h ) = E(x Bj ) / d 2 b / dzp(b, z)) / dl c T h (z + l c , oo) X; C%(A, x Bj ,Q 2 )D^(z h , l c ) (2) 

J J —oo JO 

[Fi(z h )T q (z, z + l c )+ [ Z+lc dz 'a c q (z ')p(b, z ')Fi{z h ')T q (z \z + l c )\ 

J z 

and respectively for deuterium: 

D h D (z h )= / dlJ [ Y J C h fl {2,x Bv Q 2 )D c t {z h X)F i {z h )] (3) 

J • i 
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in (2) [] and (3) the following notations are used: 



T q (h) 



a,b) = exp(- dzia q ( h) (zi)p(b,zi)) (4) 

J a 



where: b - is the impact parameter; 
(b,z) - is the coordinate of point in which the DIS takes place on one of internuclear 
nucleons; 

(h,z ) - is the point at which the color interaction of the (anti)quark located at the 
end of string with one of the internuclear nucleons being on trajectory of movement of 
knocked out quark takes place; 

E(x B j) - function, which takes into account the EMC-effect, the suppression of nuclear 
structure function at large values of the x B j P"2"| l; 
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p(r) - is nuclear density function. To parametrize this density the shell model Jl5 
used for light nuclei , and the Woods-Saxon parametrization is used [0] for middle and 
heavy nuclei. The normalization condition is Jdrp(r)=A; 

Zh =Zhv/v(t), description of z/(t) will be given after expression (6); 
n - is the ordinal number of hadron produced from the string. The hadron created at 
the fast end of string (usually it is called the leading hadron) has the ordinal number n 
= 1, hadron created behind it has the ordinal number n =2 and etc. To simplify the 
numerical calculations the limited values n < 3. were used 

C^^K^xsjiQ 2 ) - are the functions, which take into account the probability of hadron's 



Strictly speaking the expression (2) is valid in valence region xbj > 0.2, since in this region leading 
parton in string is quark. In "sea" region leading parton may be also antiquark, which color interaction 
cross section can be different from it for quark cr| ^ a c q . However for the sake of simplicity we do not 
take into account this difference. 
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h creation (with given type, charge and ordinal number i in string) from quarks and an- 
tiquarks being into the internuclear nucleon on which the DIS takes place and qq pairs 
created in color field of the string. The detailed shape for this function will be presented 
in Appendix A. 

Df(zh,l c ) - are constituent formation length l c distribution functions for hadron i from 
string which carry away the partial energy Zh- Detailed shape for these functions is given 
in Appendix B. 

Fi(zh) - are probabilities that i-th hadron produced from string carry away the part of 
the virtual photon energy Zh- Detailed shape for this function is given in Appendix C. 
o"/i(z) is the hadron-nucleon total cross section. 

cr q (z) is the total cross section of quark- nucleon interaction which could be presented 
as a sum of two components. 



a q (z) = a\{z) + a° q (z) (5) 

where (Jq(z) is part of the total cross section, which is connected with color interaction 
without essential change of energy and momentum of interacting partons. This part of 
cross section can be parametrize according to fll2|: 



z '-z- Q 2 {t) = gg^; Kt) = z/-kt; C = 1.32 mbGeV 2 . 



Parametrization (6) has bottom limit of validity which is equal Q 2 = Qo = 0.06 
Gel/ 2 . At Q 2 < Q 2 a^Q 2 ) = C/(QQ ) is using. The shape of (6) is obtained from the 
idea of color transparency, therefore naturally suppose that other part of quark nucleon 
cross section has the the same shape: 



"5W = m 

where C ' is the constant generally speaking different from C. However in order to es- 
cape the supplementary fit we assume in this work that C = C. Really in hadronization 
process we have deal not with isolate quark, but with quark which together with antiquark 
compose the string or piece of string (hadron). This system is colorless quark-antiquark 
dipole. Consequently, instead of (7) we can use the phenomenological dipole cross section. 
For example simple and convenient parametrization has been suggested in |23]. Easy to 
see that for the HERMES v and Q 2 kinematic region the expression (9) from turn 



in formula (7) of our work. Parameter C ' which obtained in this manner from dipole 
cross section is weakly changing function over v and Q 2 , with the middle value of C 
~ 1 mb-GeV 2 . However, taking into account, that in energy range W 2 = 4 50 GeV^ 2 



the value of a from p4j seems to be underestimated, we are using a = 25 mb, and 



obtained for C numerical value close to C. 

In fact to describe the attenuations measured experimentally as function of is, Q 2 , or 
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Zh, we are integrating the expression (2) ower two of three mentioned above variables. 

The calculations were performed using the following values for total cross section of 
hadron h with internuclear nucleon (we suppose that cross section for the interaction of 
hadron with proton and neutron are equal): 

a 7T + = a n - = 25 mb; = 17 mb; o^- = 23 mb. 

Some authors are using the inelastic cross section instead of total ones Unfor- 
tuntaly, the available experimental information does not allow to make a choice between 
these two posibilities. 

In present work the value of 1 GeV/fm is used for the string tension. 

Also for all calculations the following kinematical conditions were used: 

P h > 0.5 GeV, Q 2 > 1 GeV 2 , v > 4 GeV, W 2 > 4 Gel/ 2 , y < 0.85, z h > 0.1 
where Ph is the hadron momentum, W 2 is the hadron invariant mass square, y—u/Eo, Eq 
is the incident lepton energy {Eq = 27.5 GeV). 

With the presense of more data it will be possible to make a global fit to determine the 
optimal values of string tension, the parameter C in quark-nucleon cross section as well 
as to precise the value of a n ± and ak± cross-sections. 



2. Results and Discussion 



The obtained results for the nuclear attenuation of charged pions and kaons (R^ k ^) 
calculated for four different targets(N, Ne, Kr, Xe) are presented on the figures 1-6. It 
is clearly seen that in case of pions attenuation as a function of v there is no difference 
for opposite charges of pions. In case of kaons this difference becames noticeable at low 
values of v {y < 12 GeV). Also one can see that the attenuation for 7r(K)~ in case of 
symmetric targets is larger than for 7r(K) + . While in case of non symmetric targets the 
situation is vice versa. For z% dependence one can note that the difference between op- 
posite charged hadrons is increased for pions, particularly for heavy targets. Because 



of mainly HERMES Collaboration || £l|] has produced the data with different charged 
pion's attenuations during last years, the calculated curves were tuned to the HERMES 
kinematics. The following additional cuts were applied: v > 7 GeV, z h > 0.2, 4.0 < Ph 
< 13.5 and 4.0 < Ph < 15.0 GeV for N and Kr, respectively. One can see good enough 
agreement for the predicted //-dependence of n , attenuations in case of N and Kr (Fig. 
7a, c). For Zh dependence only sum of 7r + and 7r~ attenuations has been published and 
except one point there is also good agreement for calculated values and data (Fig. 7b). 
In conclusion, we have to note that developed approach can provide much better de- 
scription of data after the fitting procedure based on all existing data for different type 



of hadrons and targets. The parton distributions used for this approach |T8[ could be 
modified using the last results on parametrization for special case of nuclear parton dis- 
tributions. As it was discussed above in the valence region the quark is knocked out from 
the nucleon, so the color string is streched between this quark and nucleon remnant, and 
the final hadrons are produced from the string. In the "sea" region mainly the "sea" 
(anti) quark is knocked out from the nucleon and two strings are created. One of them 
with the knocked out (anti) quark takes away the main part of virtual photon energy 
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passing the next string an energy AE which is enough to produce 1-2 hadrons or one 
resonance. In the framework of Regge model such production of hadrons from the big 
and little strings is called "Undeveloped Pomeron". This mechanism could lead to some 
changes of the attenuation values (particularly at low u), as well as to relative changes 
for different hadrons and charges for the same hadron. As it is followed from the Regge 
model the AE value is estimated to be 1 - 3 GeV and this value also should be fitted. The 
preliminary estimations of possible influence of the mentioned above scheme with AE ~ 
1-3 GeV showed that they expected changes in obtained results will be not essential. 
In generally we have to note that the attenuation values for different types of hadrons 
as well as for the oppositely charged mesons of the same kind are very close. For the 
most heavy target 131 Xe the difference between the charged mesons is less than 5 %, and 
for different mesons is less than 10 %. In the same time one can see strong enough v 
dependence for attenuations within the HERMES kinematics range. 

Acknowledgments 

We would like to thank our colleagues from HERMES collaboration for useful dis- 
cussions. 



7 



Appendix A: Functions Cf { (A,XBj 

Notations: i = 1,2,3 means the ordinal number of three hadrons produced last in 
string(and the first number corresponds the last one) 

a) 7r + (ud) 



Z 4 1 - 

C]i = {^[-(u v (x Bjl Q 2 ) + u s (x Bj ,Q 2 )) + -d s (x Bj ,Q 2 ) 

+-^[g( d v( x Bj, Q 2 ) + d s (x B j, Q 2 )) + ^u s (x Bj , Q 2 )]}% 



where A - is atomic number; Z and N - numbers of protons and neutrons in nuclei; 
7 g - is the probability of light qq pairs (underproduction in color field fI7|| ; u v (x B j,Q 2 ), 



d v {x B j,Q 2 ), u s (x B j,Q 2 ) etc. - parton distribution functions in proton. We used parametriza- 



tion of NLO(MS) parton distributions from [|I~8|j 



Z 4 1 
C n = il[ln( u v + u s + u s ) + „{d v + d s + d s + s s + s s )) (9) 
J a y y 

+-j[7:(dv + d s + d s ) + \{u v + u s + u s + s s + s s )}}j 2 
a y y 

b) 7T-(du) 

Z 1 4 N 1 4 - 

C n ^{j[g(^ + d s)+gU s ) + ^[-(u v + U s ) + -d s )]}j q (10) 

Cj 2 and CJ 3 the same as in previous case. 

c) fc + (us) 



ZA 



N A 



Cfl = -jig( U v + U shs + gSsTg] + -j[g( d v + d s)ls + g^Tg] 



°/2 - °/3 - °/2 7? 

where 7 S - is the probability of ss pairs production in color field |T7 
d) &~(su) 



111 



C fl = j\-g Ss ^ + g Us ^ + ^[g Ss7g + g 6 ™ ^ 12 ) 



2 In further (for brevity) we omit arguments of parton distribution functions 
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°/2 _ °/3 — U /2 

For the all described here hadrons: = Cj 2 , where i = 3, 4, 5 

Appendix B: The distribution of constituent 
formation lengths in Field- Feynman fragmentation scheme 

Df(x,l c ) - is the distribution function of constituent formation lengths l c , for i-th 
hadron in string carry away parth of energy x of virtual photon. To obtain these functions 
in it was used the function f(x), which have the meaning of probability, that the first 
hierarchy (rank - 1) primary meson leaves the fraction of momentum x to the remaining 
cascade. We are using function f(x) obtained in [^0|: f(x) = 1 - a + 3ax 2 



As it is well known the it and K mesons could be produced in direct way and also as a 
result of resonances decay. In this work only the direct production is taking in account. 
In case of pions it was estimated the contribution coming from the p meson decay and 
it was shown that the second mechanism (via resonances) did not change the values of 
calculated attenuations. 

As noticed above, we are limited only with three hadrons on the fast end of string. For 
these hadrons the distribution functions are given as 

D c 1 (x,l c ) = f(l-x)5[l c -(l-x)L] (13) 

D\U.I,) = — U- /(_L)±/( x +2) (I!) 
X + t x + f L> L 



l_ ¥ , l I, 

L 

lc 



I — If 1 din 1 A- tTj 

DtM = ^f{-^)j ^f( Ul )f(^t) (15) 

X + f X + f LJ Ul min U\ LU\ 



where ui min = ' c + xL ; L - length of string (L = v/k). 
The general formula for Df(x,l c ) is given as 



L s 1 f l dUl el s f l dUi-2 ,., \ l c + xL 



Dt(x,l c ) = -—77/(^17)7 / ^/(tii).- / ^/te-2)/( /: - ) (16) 

X + -J- X ~\- ~jr Ll Ju\ min U\ Jlti-2min «i — 2 LlU\...Ui—2 

where i = 3, 4, 5, . . ., and Ui min = T lc+xL — 



Appendix C: The functions i^(x) 



Fj(x) - is the probability, that the i-th hadron in string carry away the fraction of 
momentum x. 

The general formula for i^(x) is given as 

f.W=/'^/W/' ^/W-../' x —) (it) 

Jx T]i Jx/ni V2 Jx/vu.-m^ T]i-\ VlV2---Vi-l 
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For cases of i = 1, 2, 3 

F 1 (x) = f(l-x) 
where f(x) was definded in Appendix B 

f 2W =/;^/(,)/ ( i-£) 

F3W = / 1 f/W/'fi/(, 1 )/(l 

Jx lj Jx/rj ijl 
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Figure 1: The v- dependence of the ratio R\ for the different targets: a) 14 N; b) 20 Ne; 
c) 84 Kr; d) 131 Xe. The dashed curves correspond to the case -k~ and dotted - n + mesons 
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Figure 2: The same as in Fig. 1, but calculated for the charged kaons 
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Figure 3: The - dependence of the ratio R\ for the different targets: a) 14 N; b) 20 Ne; 
c) 84 Kr; d) 131 Xe. The dashed curves correspond to the case 7r~ and dotted - n + mesons 
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Figure 4: The same as in Fig. 3, but calculated for the charged kaons 
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Figure 5: The Q 2 - dependence of the ratio R\ for the different targets: a) 14 N; b) 20 Ne; 
c) 84 Kr; d) 131 Xe. The dashed curves correspond to the case 7T~ and dotted - tt + mesons 
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Figure 6: The same as in Fig. 5, but calculated for the charged kaons 
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Figure 7: The comparison of predicted values for attenuation with HERMES data: a)^ 
vs v for 14 iV b) 7i + + 7r~ vs Zh for 14 iV j^; c)tt vs v for S5 Kr ffi^ (preliminary results) 
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